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Abstract. In this paper we present the approach we have used to determine and account for the anisoplanatic point 
spread function (PSF) in deep adaptive optics (AO) images for the Survey of a Wide Area with NACO (SWAN) at 
the ESO VLT. The survey comprises adaptive optics observations in the K s band totaling ~ 30 arcmin 2 , assembled 
from 42 discrete fields centered on different bright stars suitable for AO guiding. We develop a parametric model 
of the PSF variations across the field of view in order to build an accurate model PSF for every galaxy detected in 
each of the fields. We show that this approach is particularly convenient, as it uses only easily available data and 
makes no uncertain assumptions about the stability of the isoplanatic angle during any given night. The model was 
tested using simulated galaxy profiles to check its performance in terms of recovering the correct morphological 
parameters; we find that the results are reliable up to K a ~ 20.5 (Kab ~ 22.3) in a typical SWAN field. Finally, 
the model obtained was used to derive the first results from five SWAN fields, and to obtain the AO morphology 
of 55 galaxies brighter than K s — 20. These preliminary results demonstrate the unique power of AO observations 
to derive the details of faint galaxy morphologies and to study galaxy evolution. 
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1. Introduction 

In recent years an increasing number of adaptive optics 
(AO) systems have become available, which are capable 
of obtaining near-infrared images at the diffraction lim- 
its of 8-meter class telescopes. The scientific potential of 
such instrumentation is evident, but careful analysis of the 
resulting data is necessary to take account of the effects 
introduced by the method and the limitations of wavefront 
sensing and correction. A key point in the detailed analy- 
sis of AO images is the determination of the point spread 
function (PSF) across the whole field of view. However, 
this task is made more difficult because the performance 
and correction of the AO system are strongly anisopla- 
natic, and depend on many parameters such as the bright- 
ness of the reference guide star and the structure (height, 
size-scale, velocity) of the often quickly changing atmo- 
spheric turbulence. As a result, the PSF produced by an 
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AO system can change rapidly in both time and position 
on the frame. For natural guide star (NGS) wavefront 
sensing, an extragalactic science target will typically be 
offset from the guide star by A9 ~ 10"-30" at best. As a 
result, the on-axis PSF no longer provides a suitable refer- 
ence model for the off-axis PSF at the position of the sci- 
entific target. Until the advent of multi-conjugate adaptive 
optics systems on 8-m class telescopes (e.g., Le Louarn et 
al. l2002fl full correction over a wide (1-2') field will not be 
possible, and even then there will still be some PSF varia- 
tion across the field (e.g., Verinaud et al. 2003). Therefore, 
a simple and generally applicable technique to model the 
variations in the PSF across a given science field would be 
of benefit to all wide-field adaptive optics data. 

In this paper we present the approach we have used 
to account for the anisoplanatic PSF in deep adaptive op- 
tics images for the Survey of a Wide Area with NACO 
(SWAN). In the following section SWAN will be briefly 
described and motivated. In Sect. [3] some different pro- 
posed methods for PSF reconstruction will be analyzed, 
while in Sect. 0] our method will be explained and dis- 
cussed. Computer simulations were used to validate the 
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model, and their results are shown in Sect. [S] The first 
application of the model to five of the deep science fields 
will be presented in Sect. [HI and our conclusions follow in 
Sect.d 

2. SWAN 

In order to make full use of the new generation of tele- 
scopes, it is necessary to overcome the blurring effects of 
the atmosphere through the use of AO systems. These 
can allow ground-based telescopes to operate at or near 
the diffraction limit in the near infrared (~ 0.06" in K- 
band for an 8 meters telescope), resulting in high angular 
resolution and a low background in each pixel. In princi- 
ple, such capabilities should offer important benefits for 
studying how galaxies form and evolve in the early uni- 
verse. In practice, however, reaping the expected rewards 
has proved difficult, due to the small number of extra- 
galactic sources known lying at distances Ad 30" from 
bright (V Si 13) stars suitable for AO guiding. 

The prospects for AO cosmology will undoubtedly im- 
prove with the widespread adoption of laser guide star 
(LGS) systems, since these impose less stringent require- 
ments on the brightness and nearness of stars used for 
tip-tilt correction. In the mean time, to overcome the cur- 
rent shortage of extragalactic AO targets for NGS wave- 
front sensing, we have undertaken a program to identify 
and characterize faint field galaxies lying close to bright, 
blue stars at high Galactic latitudes (see, e.g., Larkin et al. 
1999] Davies et al. l2001J) . Our own 42 southern bright star 
fields were initially imaged at seeing-limited resolution 
in K s with SOFI at the ESO New Technology Telescope 
(Baker et al. l2003[) . were followed up with optical imaging 
(Davies et al. 2005), and are now targets for VIMOS opti- 
cal spectroscopy at the ESO Very Large Telescope (VLT). 

SWAN is the AO-assisted payoff of these seeing-limited 
preliminaries. Having already characterized large samples 
of objects in our bright star fields, we targeted them with 
NACO on the VLT in order to exploit the present gen- 
eration of AO technology for galaxy evolution studies. 
NACO comprises the NAOS Shack-Hartmann AO mod- 
ule (Rousset et al. 12003)) mated with the CONICA near- 
infrared camera (Lenzen et al. H998fl . Our choice of NACO 
observing mode was dictated by our desire to differentiate 
SWAN from previous HST/NICMOS surveys. First, we 
chose to prioritize survey area over depth, thereby improv- 
ing SWAN's sensitivity to rare objects and its robustness 
against cosmic variance. Second, we chose to image in K s , 
where NICMOS is less sensitive than in J and H, thus 
making SWAN preferentially sensitive to red objects. Use 
of NACO's 0.054" pixel scale (to maximize field of view) 
and the Strehl ratios of 30-60% typically achieved in K s 
thus result in images that are undersampled. Each NACO 
pointing provides only a usable ~ 0.75 arcmin 2 of the full 
55.5" x 55.5" detector area, due to losses from dithering 
and the central star (see, e.g., Fig. Nevertheless, the 
anticipated survey area that will result from assembling 




Fig. 1. Example of a SWAN field: SBSF 18. The circles are 
the extended objects detected by SExtractor (SExtractor 
stcllarity index SSI < 0.9), while the squares are point 
sources (SSI > 0.9). A ghost of the bright guide star is 
marked with a cross. 



42 such images will be - at ~ 30 arcmin 2 - some six times 
larger than the NICMOS survey of the HDF and flanking 
fields in J and H (Dickinson Dickinson et al. 15000)1 . 

In the ecology of near-IR surveys, SWAN aims to occupy 
a niche combining the high angular resolution of a space- 
based survey with the shallower depth and wider area of 
a ground-based survey, thereby probing sources that are 
compact, faint, red, and rare more effectively than any 
other survey to date. First results from an initial analy- 
sis of nine SWAN AO fields are presented in Baker et al. 
( 2005), where NACO imaging is seen to detect a popula- 
tion of compact galaxies that cannot be identified as such 
in seeing-limited data. 

To extract full information about the morphology of 
the sources detected in the SWAN images by taking into 
account the effects introduced by the PSF, it is necessary 
to develop a reliable model of the off-axis PSF in each of 
the fields. This task is not easy, however, as the AO PSF 
changes quickly in both time and position on the frame; in 
our case it is made even more difficult by the distinctive at- 
tributes of the SWAN observing strategy. First, because of 
their high Galactic latitudes, relatively few point sources 
are present in each of the fields, so that very few objects 
can be used as references to constrain how the PSF varies 
off-axis. In addition, in order to be background-limited so 
that the aim of the survey to detect faint galaxies can be 
realized, the on-axis star is always saturated on the science 
frames, meaning the on-axis PSF has to be obtained with 
separate unsaturated exposures. Finally, NACO has only 
very rarely been used to image stellar fields with the (un- 
dersampled) 0.054" pixel scale and the visible wavefront 
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sensor; as a result, the VLT archive contains no suitable 
stellar fields that can help characterize the properties of 
the PSF as far off-axis as the SWAN images extend. 

Given the above considerations, our data obviously 
pose several challenges for the recovery of intrinsic source 
parameters. Indeed, the same challenges will be faced by 
any extragalactic AO survey that relies on NGS or LGS 
wavefront sensing, and that efficiently builds up statisti- 
cal samples of faint field galaxies by including a wide area 
around each bright star position. It is thus of general in- 
terest to devise a strategy to account for anisoplanaticism 
in deep, wide-field AO imaging. 

3. Other proposed methods to quantify 
anisoplanaticism 

There have been several methods proposed in the litera- 
ture for estimating or calculating the off-axis PSF in an 
AO field. In this section we summarize these methods, 
and explain why they are not ideal for deep wide-field AO 
imaging with current instrumentation. 

A theoretical analytical expression to model PSF vari- 
ation and anisoplanaticism has been derived by Fusco et 
al. H2000f) . who validated their method on both simulated 
and experimental data, leading to a reduction in the error 
on the magnitude estimation in stellar fields from more 
than 30% to only 1%. They showed that the total optical 
transfer function (OTF) is simply the product of the on- 
axis OTF with an anisoplanatic OTF. The on-axis OTF 
can be readily obtained, for example from real-time data 
accumulated by the AO system or a measurement of the 
on-axis PSF. The anisoplanatic OTF, however, requires 
an independent measurement of how the turbulence in the 
atmosphere depends on altitude, as quantified by the at- 
mospheric refractive index structure constant C„. In their 
example case this was obtained by balloon probes, but 
similar measurements are not usually available. 

In a similar vein, the possibility of using AO measure- 
ments together with simultaneous scintillation detection 
and ranging (SCIDAR) measurements to reconstruct the 

profile and hence the off-axis PSF, was investigated 
by Weifiet al. I|2()()2al l2002b)> . While such measurements 
were made and analyzed, the idea was never brought to 
fruition and no such system has been installed as a per- 
manent facility instrument. 

An alternative semi-empirical approach was proposed 
by Steinbring et al. l|2002|l . based on calibration images 
of dense stellar fields to determine the change in PSF 
with field position. Their results showed that this simple 
method reduces the error in the prediction of the FWHM 
of the PSF at large distances off-axis from 60% to only 
~ 20%. However, obtaining a suitable calibration field is 
not an easy task as atmospheric conditions can change 
rapidly. Steinbring et al. were able to complete the obser- 
vations necessary to construct their mosaic images in less 
than 10 minutes. However, this approach is not practicable 
for even moderately deep- field observations. Our SWAN 
observations, for example, typically have integration time 



(without overheads) of 60 minutes. The additional time 
needed to switch continually between science and cali- 
bration fields would be prohibitive - even assuming that 
a suitable stellar field with a guide star of comparable 
brightness can be found not too far from the science field 
and also at a similar airmass. Furthermore, the authors 
point out that variations of up to 50% in the measured 
on-axis Strehl ratio ultimately limit the accuracy of the 
results. 

A third approach to the problem was investigated by 
Tristram & Prieto l|2005[) . They parameterized the PSF 
using an elliptical Gaussian and a decaying elliptical ex- 
ponential function. Starting with 11 parameters, they were 
able to eliminate all but four, which would need to be de- 
rived from stars in the science field. At high Galactic lat- 
itudes, the difficulty is that the number of point sources 
will be very limited in any science fields. To constrain 4 
parameters reliably, the authors used 20 point sources, 
while our SWAN images always have fewer point sources 
per field - in a number of cases there are fewer than four 
- that are usually too faint for a detailed PSF fitting as is 
required here. 

Given that it is not practical to apply any of these 
methods to SWAN or similar data, in the next Section we 
propose an alternative method suited to cases where there 
are very few unresolved sources in the science field (and 
those which do exist are often rather faint) , and where one 
does not necessarily need a highly precise estimate of the 
PSF. 

4. The PSF model for SWAN 

To tackle the issue of anisoplanaticism in the SWAN data 
and to model an accurate PSF for every galaxy we de- 
tect we have developed a parametric model of the vari- 
ations of the PSF across the field of view. This kind of 
approach is particularly convenient, as it uses only eas- 
ily available data and makes no uncertain assumptions 
about the stability of the isoplanatic angle during any 
given night. Although the PSF reconstruction is not per- 
fect, we will show that it is accurate enough for our pur- 
pose of measuring galaxy morphologies. 

According to a theoretical analysis of the anisoplanatic 
effect in AO systems (e.g., Fusco et al. 2000 Voitsekhovich 
& Bara 1999 : ), the off-axis PSF can be expressed as a con- 
volution between the on-axis PSF of the guide star with 
a spatially variable kernel, i.e., 

PSF(f,0) = PSF(f,O)<8>K(f,0) (1) 

We will make the extreme assumption that the model of 
the kernel K(f , a) is sufficiently well constrained that the 
difference in the kernel from one science field to the next 
is reduced to a single parameter, namely the isoplanatic 
angle. This is the most convenient solution for fields with 
sparse point source populations and, as we show, is able to 
reproduce a reasonable approximation to the PSF at any 
given location in the field of view using just the the on- 
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Field 


Date 


Filter 


Exp. Time 


Strehl 








(min) 


(%) 


NGC 6809 


15.12.02 


K s 


60 


28 


NGC 6752 


06.09.03 




12 


27 



Table 1. Observations of PSF calibration fields. The 
Strehl ratio is that measured on-axis in the short exposure 
frames taken before and after the deep science exposures. 



axis PSF and a few reference point sources in the science 
field as calibration. 

For each SWAN field the on-axis PSF was monitored 
by taking short unsaturated images of the bright star both 
before and after (and in some cases half-way through) 
the deep science exposure. These short exposure frames 
were taken in the K s band, IB_2.21 (A c = 2.210 /an, 
FWHM = 0.060 urn) or NB.2.17 (A c = 2.166 fim, FWHM 
= 0.023 (im) filters. The integration time varied in the 
range 0.5-2.0 sec depending on the filter used. In addition, 
observations of galactic star cluster fields were planned pe- 
riodically in order to build up a database of anisoplanatic 
PSFs. At the time of writing only two such calibration 
fields had been observed (see Table 0}, but between them 
and the 7 SWAN fields with the highest numbers of point 
sources (see end of this section), included in Fig. 2, they 
still provide nearly 80 different point sources. These fields 
will be used in the following analysis as a test bench to 
derive a suitable parametric model for the off-axis PSF. 

The images were reduced using PC-IRAF version 
2.11.3. The presence of the bright star in the center of 
a field less than 1' across made the data reduction a little 
more complex then usual, requiring two iterations of ob- 
ject masking and sky subtraction to avoid over-subtraction 
from very extended faint scattering around bright objects. 
For further details about data reduction see Cresci et al. 



The isoplanatic angle of these fields was measured by 
considering - independently for each field - the Strehl ra- 
tio of the detected stars as a function of radius from the 
guide star. To these data we fit the theoretically expected 
function (e.g., Beckers ITMfl Roddier lT5§5j> : 



SR tx exp 



A9 





(2) 



where A9 is the offset angle from the guide star on the 
sky and 9q the isoplanatic angle, defined as the separa- 
tion angle at which the Strehl ratio has degraded by a 
factor e with respect to the on-axis value. Although in a 
strict sense this refers to a telescope with infinite aperture 
(e.g., Beckers Il993fl . for our implementation this defini- 
tion is adequate. A constant term was added to Eq. [21 in 
order to match the observed data. The resulting best fit- 
ting values for 6 are 9.30" and 6.61" for NGC 6809 and 
NGC 6752 respectively. These are surprisingly small given 
the mean isoplanatic angle of the other seven SWAN fields 



- 0.2 



• NGC 6752 
■ NGC 6B09 

* SWAN fields 




12 3 4 5 

Radial distance (9 ) 

Fig. 2. Isoplanatic angle scaling. The Strehl ratios of point 
sources in the two calibration fields and 7 SWAN fields 
(SBSF 14, SBSF 15, SBSF 18, SBSF 24, SBSF 27, SBSF 28, 
SBSF41) are plotted against their radial distance a from 
the AO guide star (in units of the corresponding isopla- 
natic angle). The black solid line is the theoretical curve 
used for fitting. 



of ~17". However, they do provide a larger range in 9q over 
which our model can be constrained. The radial distances 
of all the point sources in the two calibration fields and in 
the 7 observed SWAN fields could then be normalized to 
their isoplanatic angles. This is shown in Fig. [21 where the 
Strehl ratios of the point sources in these fields are plotted 
against their rescaled separations from the AO guide star. 
It can be seen that all the points now lie on the same the- 
oretical distribution, showing that the isoplanatic angle is 
the source of the principal variation between the different 
fields. 

Our model is motivated by the fact that the wave- 
front error in the anisoplanatic kernel is dominated by 
the tip-tilt terms. This is borne out in practice, since the 
major degradation of the PSF observed in the calibra- 
tion field images is an increasing radial elongation towards 
the guide star as one moves further off-axis. This suggests 
that, as a first approximation, one can represent the aniso- 
planatic kernel K(f , a) with an elliptical Gaussian elon- 
gated towards the guide star. For each of the point source 
in both of the calibration fields, we therefore derived the 
best fitting radial and tangential FWHMs of the elliptical 
Gaussian kernel to be convolved with the on-axis PSF in 
order to reproduce the observed off-axis PSF (see Fig. [3]). 
The fit was obtained by minimizing the sum of the squared 
difference between the observed star and models with dif- 
ferent kernel FWHMs, weighted by the flux in each pixel 
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Fig. 3. Upper panels: FWHM fits to the kernel. The radial (left panel) and tangential (right panel) FWHMs of the 
kernel to be convolved with the on-axis PSF to reproduce the observed PSF for the point sources in the calibration 
fields (diamonds) and in the seven SWAN fields of Fig. |2] (pluses), along with the best fitting second order polynomial. 
Lower panels: The total FWHMs of the observed point sources (radial at left, tangential at right) as a function of 
radial distance a from the guide star in units of the corresponding isoplanatic angle 9q, compared to the best fitting 
model for the kernel convolved with the guide star's FWHM. 



of the star frame (in order to optimize the fit in the PSF 
core) and including all pixels brighter than 1% of the star 
peak. As the noise is dominated by the background, which 
is constant over each star, a noise term was not included 
in the expression, because it would have added just a con- 
stant scaling factor. 

The residuals of the fitting are shown in Fig. 0] as a 
function of separation from the guide star. The residuals 
are defined as follows: 

Residuals = ( (Star, - Model )^ • S^arQ 
(EiCStar,)) 1 - 5 

where Star is extracted from the NACO data and Model is 
the best fitting kernel convolved with the on-axis PSF, so 
that the numerator is the quantity minimized during the 
fit, which was then normalized using the star brightness. 
It is clear that not all of the stars are equally well repre- 
sented by the model, as expected due to the very simple 
kernel adopted to reproduce the PSF. However, the resid- 
uals do not depend strongly on the field in which the star 
lies (see Fig. suggesting that the approach is robust to 
different conditions. The contribution of the noise to the 
residuals could be significant in the fainter point sources 
in the SWAN fields. 

The FWHMs obtained along the radial and tangential 
directions were then fitted empirically as functions of the 
radial distance from the AO guide star. The function cho- 
sen was a second order polynomial, which was forced to 



0.8 I 



0.6 - 



"w 



0.2 - 




12 3 4 5 

Radial Distance (8 Q ) 

Fig. 4. Residuals from the PSF fitting with the elliptical 
kernel for the point sources in the two calibration fields, 
NGC 6752 (triangles) and NGC 6809 (squares), and in the 
7 SWAN fields of Fig. H (stars). The residuals plotted are 
defined in Eq. [3J 



a constant after reaching a maximum. Using higher order 
polynomials (or other analytical functions) just increased 
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Fig. 5. Comparison between the true and model PSFs for 
a point source in NGC 6752. The left panel shows the 
original star, the central panel the model for the corre- 
sponding position on the detector, and the right panel the 
residuals after subtracting the two. The star is located at 
26.8" from the AO guide star. 



the number of free parameters without improving the fit. 
The final best fitting functions obtained are 
for the radial axis: 

0.107 • a - 0.00594 -a 2 if a < 9.1 

(4) 



FWHM r (") 



0.481 



if a > 9.1 



and for the tangential axis: 

0.0871 • a - 0.00951 -a 2 if a < 4.6 



FWHM t (") 



0.199 



if a > 4.6 



(5) 



where a = 9/9 is the radial distance from the guide star 
rescaled using the isoplanatic angle. We notice that, al- 
though the method could be in principle applied to every 
AO system, the derived parameters were derived specif- 
ically for our SWAN dataset, and they have to be re- 
calibrated in order to use the same PSF model for any 
different adaptive optics system, instrument setting and 
operating conditions. 

The variation of the total FWHM - i.e., the kernel con- 
volved with the on-axis PSF - along both radial and tan- 
gential directions for the point sources in the calibration 
field is shown in the lower panels of Fig. EI It can be seen 
how the total FWHM has a behavior consistent with what 
is observed for other AO PSFs (e.g., Flicker & Rigaut 

HE2Ji. 

Using the derived parametric model for the kernel, it 
is now possible to build the model PSF for any position 
in any of the SWAN fields with only an estimate of the 
isoplanatic angle in the respective field. The isoplanatic 
angle is derived by fitting the Strehl ratio of the few point 
sources detected in each science frame as a function of 
radius with Eq. We have used SExtractor (Bertin & 
ArnoutsHnMIl to identify the sources in our SWAN NACO 
images, and as reference point sources for fitting the Strehl 
ratio, we use those with SExtractor stellarity index SSI > 
0.9 (see Fig0. The derived isoplanatic angle is then used 
to rescale the radial distance from the guide star and to 
build the model PSF. 

We additionally compared these results to those that 
could be obtained using PSFs generated by the PAOLA 
software (http:/ /cfao. ucolick.org/software/paola.php), with 
an appropriate isoplanatic angle. PAOLA (Performance 
of Adaptive Optics for Large Apertures) is a set of func- 
tions and procedures written in IDL for calculating the 



performance of an AO system installed at the focus of 
an astronomical telescope. It relies on an analytic expres- 
sion of the power spectrum of the corrected phase and its 
relation to the AO OTF. It therefore assumes that the 
AO system is "perfect", failing to reproduce second or- 
der effects due to imperfections of the system. A set of 
PSFs were generated at different distances from the guide 
star, using in the model average parameters for the VLT 
and NACO. The radial distances of the PSFs were then 
rescaled to match the isoplanatic angle observed in the dif- 
ferent fields. However, the match with the true PSFs was 
no better - and often worse - than that obtained with our 
parametric model, in particular worsening for increasing 
distances from the guide star. 

5. Testing the PSF model with simulations of 
galaxy profiles 

The PSF model developed in the previous section was 
carefully tested to see how well it performs in terms of 
recovering the correct morphological parameters of sim- 
ulated galaxy profiles. The galaxy models were built by 
convolving a Sersic l|196 8) profile, 



J(i?)-/(i? e )xe(- fc "X[^) 1/ "- 1 ]) 



(6) 



(where R e is the effective radius, n is the Sersic index and 
b n is a constant that varies with n) with a true PSF ex- 
tracted from either a calibration or a SWAN field. Our 
tests are particularly robust since we can compare the pa- 
rameters recovered using the model PSF not only with 
those used as input for the simulations, but also with those 
recovered using the true PSF used to convolve the simu- 
lated profile. The primary aims of the simulations were to 
determine how well it was possible to estimate as a func- 
tion of magnitude (1) the Sersic index n - and in particular 
whether it was possible to separate disk and elliptical pro- 
files, and (2) the effective radius R e of the profile, which 
indicates the size scale of the galaxy. 

We used Sersic index n = 1 for disk-like galaxies and 
n = 4 for elliptical-like galaxies. For both types we ran 
several sets of simulations, each including 100 galaxies, at 
fixed magnitudes ranging from 17-22 mag in the K s band, 
with the same pixel scale used in the SWAN images. The 
noise was set to the level expected for NACO integration 
time of 60min for all the galaxies, since this is the typical 
integration time for a SWAN field. The inclination and po- 
sition angle were random for each simulated profile, while 
the distribution of effective radii R e was chosen to roughly 
reproduce that observed in the real data - i.e. 40% of the 
objects having R e = 0.1", 20% each with R e = 0.2" and 
R e = 0.3", and 10% each with R e = 0.4" and R e = 0.5". 
Ten different PSFs extracted from point sources in the 
calibration or SWAN fields were used for each set of 100 
galaxies. 

To recover the galaxy morphological parameters, we 
used GALFIT (Peng et al. I2002[> . a widely used software 
package which fits an image of a galaxy and/or point 
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Disks 



Ellipticals 





Mag 

Fig. 6. The left panels show the RMS of the magnitude, Sersic index and effective radius obtained from the GALFIT 
fits of the simulated profiles with Sersic index n — 1 as a function of the input magnitude. For each point, 100 galaxies 
with 10 different PSF were used (see text for details). The solid lines are the RMS obtained using the model PSF as 
an input for GALFIT, while the dashed ones are obtained with the original NACO PSF. The right panels show the 
results obtained for simulated profiles with Sersic index n = 4. In both cases it is clear that the parameters can be 
reliably recovered up to K s ~ 20. 



source with one or more analytic functions. For our sim- 
ulations we used GALFIT to fit a single Sersic profile, 
leaving as free parameters the center of the galaxy, the 
Sersic index n, the effective radius, the magnitude, the 
position angle and the axis ratio. The software needs an 
image of the PSF as input in order to fit the galaxy pro- 
file, both to deconvolve the original galaxy image and to 
convolve the derived model. For each simulated galaxy we 
ran GALFIT using both the true PSF from the NACO 
data - i.e., that which had been used to create the profile 
- and the corresponding model PSF, in order to compare 
the results. 

The results of the simulations are shown in Fig. H3 
where the uncertainties in the derived magnitude, Sersic 
index, and effective radius are plotted as a function of 
magnitude for both disk-like and elliptical-like profiles. 
The results obtained using both the model PSF (solid 
lines) and true PSF (dashed lines) as input to GALFIT 
are shown, and it can be seen that these are in general 
very similar. Since using the true PSF is the best that 
can be done, we conclude that - in the application here of 
deriving the morphological parameters of galaxies - our 
model is a reasonable and valid approximation to the true 
PSF. 

The most significant discrepancy between the two sets 
of results is a larger RMS in the Sersic index for the sim- 
ulated disk galaxies when the model PSF was used in 
GALFIT. In addition, the RMS for the Sersic index is 



quite large (~ 1) for elliptical galaxies, but this time for 
both the model and true PSFs. In principle this rather 
large uncertainty on the Sersic index could restrict our 
ability to separate the two populations of disk-like and 
elliptical-like galaxies. However, in practice, the large vari- 
ance for the ellipticals is due to a long tail to high values 
of n, as apparent in Fig. [7| We have found that starting 
with two populations of simulated galaxies with Sersic in- 
dex n = 1 and n = 4 respectively, only a small fraction 
(that increases with the magnitude of the galaxies, see 
Fig. |SJ of disks was recovered with n > 2, and few ellip- 
ticals were recovered with n < 2. The two classes can be 
clearly distinguished up to K s — 20.5, where 90% of the 
disks are recovered with Sersic index n < 2 and 90% of 
the ellipticals are recovered with Sersic index n > 2, using 
the model PSF in GALFIT (see Fig. 01. As a result, it is 
possible to set a threshold of n = 2 that is able to dis- 
criminate between the two populations with a high degree 
of success. The fraction of disks recovered with a Sersic 
index n < 2 and the fraction of ellipticals recovered with 
n > 2 are shown in Fig.[S]as a function of the magnitude, 
using both the model and true PSF. As the correct fitting 
of the Sersic index is limited by the signal to noise of the 
sources, in order to correctly discriminate between disks 
and ellipticals at a fainter magnitude K s , the integration 
time should be scaled as At (hr) = io[°- 8 (^- 20 - 5 )] . 

As an additional check of the usefulness of the model 
PSF, we repeated the experiment described above but giv- 



8 



G. Cresci et al.: Accounting for the anisoplanatic PSF in deep AO images 



Mag(K)=20 



100 

50 






r-l- 


1 ' 


— 1— 


1 . 






Disks 




Ellipticals ~ 


1 









1,17 



1 2 3 4 5 

n Sersic (Input distribution) 




2 3 4 5 

n Sersic (from True PSF) 



40 
20 



1 --H 1 

~ i 

'. i 1 
_ i Disks i 


i | i | i | i | _| 
Ellipticals ^ 







1 2 3 4 5 6 

n Sersic (from Model PSF) 

Fig. 7. Disk/elliptical separation using GALFIT measure- 
ment of Sersic index for simulated profiles with K s = 20. 
The upper panel shows the input distributions of Sersic 
index, while the central panel shows the results using the 
true NACO PSF used to build the profile as an input for 
GALFIT, and the lower panel what is obtained with the 
model PSF. The discriminating value for the Sersic index 
n = 2 (see text) is shown as a gray line. 




17 IB 19 20 21 
Mag 



IS 19 20 21 
Mag 



Fig. 8. The fraction of ellipticals recovered with n > 2 (left 
panel) and the fraction of disks recovered with a Sersic in- 
dex n < 2 (right panel) as a function of the input magni- 
tude. The dashed lines show the results obtained using the 
true NACO PSF, while the solid lines the ones obtained 
with the model PSF. 



Field 


Date 


Filter 


Exp. Time 


CE 


do 








(min) 


{07 \ 
(%) 


( ) 


SBSF 14 


15.12.02 


K s 


60 


62 


11.9 


SBSF 15 


17.12.02 


K s 


60 


66 


21.8 


SBSF 18 


17.12.02 


K s 


40 


68 


21.2 


SBSF 24 


21.03.03 


K s 


60 


36 


11.8 


SBSF 41 


14.06.03 


K s 


60 


53 


19.7 



Table 2. Observations of the five SWAN fields included in 
the analysis here. CE is the encircled energy reported by 
the adaptive optics system in real time. The isoplanatic 
angle 8 is that derived by the fit using point sources in 
the fields (see Sect.0J. 



ing GALFIT the observed on-axis PSF instead of either 
the model or true PSF. The outcome was that, even at rel- 
atively bright magnitudes K s ~ 17, all the galaxies with a 
Sersic index n = 4 were actually recovered with an index 
of n < 2, making it impossible to discriminate between 
the two populations. This emphasizes the important re- 
sult of our simulations that it is only possible to discrimi- 
nate between disk galaxies and elliptical galaxies using an 
anisoplanatic PSF, and that using our very simple PSF 
model is nearly as good as using the true PSF. 

6. First results from five deep science fields 

Prompted by the results on the simulated galaxy pro- 
files, we have begun using our PSF model coupled with 
GALFIT to derive morphological parameters of the galax- 
ies detected in a few of the SWAN survey fields in or- 
der to gain some insight into what we should expect, al- 
though clearly a full analysis requires the results of the 
dedicated redshift survey currently underway. In this pa- 
per we present the first results obtained from the first five 
fields, SBSF 14, SBSF 15, SBSF 18, SBSF 24 and SBSF 41, 
all of them already included in the preliminary discussion 
in Baker et al. l|2005jl . The observations are summarized 
in Table E3 and the data were reduced using the same 
procedure adopted for the calibration fields (see Sect.0J). 

The performance obtained in these first SWAN fields 
compares favorably with that obtained by previous AO ob- 
servations of faint field galaxies, e.g., Larkin et al. |2000) 
and Glassman, Larkin, & Lafreniere (2002). They ob- 
tained 12 ff-band (1.6 /im) images of disk galaxies at 
z ~ 0.5 using the AO system on the Keck II telescope. 
They obtain FWHMs ranging from 0.050" to 0.15", and 
Strehl ratios of l%-20%. Due to the very small field of 
view (4.5") only one galaxy could be observed for each 
frame. The PSF was reconstructed using dedicated obser- 
vations of faint off-axis stars with comparable offset, but 
different position angle, from the guide star. Due to the 
very low S/N obtained, the extraction of morphological 
parameters was run under several assumptions to limit the 
number of free parameters in the fit of a bulge and a disk 
to each galaxy. In addition, the AO galaxy profiles only 
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had reasonable S/N out to 1"~2" radius, thus limiting the 
constraints on the outer portions of the disks. Steinbring et 
al. (2004) imaged with the same AO system three galaxies 
selected because they had been observed previously with 
Hubble Space Telescope at optical wavelengths. However 
they conclude that higher S/N was needed in the AO im- 
ages in order to match the HST data, and the results were 
limited due to poor constraints on the AO PSF. 

6.1. Source detection and modeling 

In each SWAN field, sources were detected using 
SExtractor, with the appropriate parameters set to pro- 
vide a positive detection for objects brighter than 1.5 a 
per pixel over an area of more then 3 pixels. To improve 
the detection of faint sources we used a Gaussian filter 
(<7 = 1.5 pixels) to smooth the image. False detections 
at the noisy borders of the mosaic and on the spikes and 
the ghost of the bright guide star were removed. For the 
former, a mask which indicted the fraction of the total 
integration time spent on each pixel was used; objects de- 
tected in pixels below a specified threshold were rejected. 
For the latter appropriate object masks were created. The 
resulting coverage of these 5 fields was 3.2 arcmin 2 , within 
which a total of 178 sources were detected down to a mag- 
nitude of K s ~ 23 {Kab ~ 24.8). Of these, 146 were al- 
located a stellarity index SSI < 0.9 by SExtractor, and 
hence classified as galaxies. The SExtractor classifications 
should be treated with caution since they assume a con- 
stant PSF across each field. However, all the objects clas- 
sified as stars by SExtractor lie on an upper envelope in a 
Strehl versus radial distance plot, i.e., they have the high- 
est Strehl ratio among the sources at the same distance 
from the guide star, supporting their classification as point 
sources. 

We again used GALFIT to derive the morphologi- 
cal parameters of the detected galaxies. We used sin- 
gle component Sersic profiles, providing the output from 
SExtractor as the initial guesses for magnitude, position, 
position angle, and axis ratio. For the deconvolution, a 
model PSF was created using the position of the center 
according to SExtractor and scaled according to the iso- 
planatic angle, which is shown for each field in Table H3 
For each object, GALFIT was run twice, using as initial 
guesses for the second iteration the outputs of the first 
iteration. 

Some of the sources still showed bumps in the center or 
disk-like structures in the residuals after the single com- 
ponent fit. These objects were re-fitted using either two 
Sersic components (bulge + disk) or a Sersic component 
with an additional point source, to improve the fit and 
minimize the residuals. 

A second iteration was also performed for close pairs or 
interacting galaxies; we fit both of them simultaneously to 
avoid incorrect background estimation or contamination 
from the companion. Some examples of single and multiple 
component fits are shown in Fig. For galaxies fitted 
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Fig. 10. Distribution of effective radii R e for the 55 
sources detected down to K s = 20 in the five SWAN fields 
observed with NACO (black hatched histograms), and for 
the 92 sources down to K s = 20 observed in in larger fields 
centered on the same stars with SOFI (gray shaded his- 
tograms). The distribution is not corrected for complete- 
ness. For galaxies with more than one component, only 
the parameters of the larger scale (i.e., non-nuclear) one 
are used. 

using multiple components, we report only the values for 
the larger scale component; the nuclear properties will be 
discussed elsewhere. 

6.2. Source morphologies 

The results obtained for the effective radii of the 55 galax- 
ies brighter than K s — 20 are shown in Fig. QJJI While 
the redshifts of these objects are presently unknown, the 
magnitude-redshift relation of Cowie et al. (1996) and the 
K20 survey (Cimatti et al. l2H0l)l indicate that at K = 20 
the median redshift is z ~ 0.8 — 1. At this redshift, our 
spatial resolution of 0.1", which also corresponds to the 
smallest effective radius bin, is equivalent to only 500 pc 
for typical cosmologies, hinting at the exciting potential 
of this work. 

Using n — 2 to discriminate between disks and ellip- 
ticals, we find there are 24 elliptical-like and 31 disk-like 
galaxies, with uncertainties below 10% on both numbers 
(see SectionE}. In fact, with similar numbers of each type 
found, the number of ellipticals scattered into the disk-like 
category should roughly compensate the number of disks 
scattered out of it, and vice- versa. 

The distribution of effective radii is strongly peaked to- 
wards very compact sources, even down to scales of 0.1". 
To check if this is due to a selection effect introduced by 
the AO correction - since adaptive optics data are inher- 
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Original Single component fit Two component fit 



Fig. 9. Examples of single and multiple components fits by GALFIT of SWAN galaxies using the model PSF. In all 
images North is up and East is left. The distance of the galaxy from the guide star is shown in the upper right corner 
of each panel. For both galaxies, the K s magnitude and the source identifier according to Baker et al. (20(31) are 
reported. Upper panel: single component Sersic fit (K s = 19.2; R e = 0.4"; n = 0.6; b/a = 0.6; PA = 75). The left 
image is the original galaxy, the GALFIT best fitting profile is in the middle, and the residuals of the subtraction 
are on the right. Lower panel: comparison of the residuals between a single component fit (in the middle) and a two 
component fit (on the right), while the original galaxy is shown on the left. The final best fit was obtained with a disk 
component (K s = 19.6; R e = 0.4"; n = 1.0; b/a — 0.7; PA = 76) and a point source (K s = 21.4) in the center. 



ently most sensitive to compact sources - we ran GALFIT 
on the same five fields, but using the seeing limited obser- 
vations made with SOFI (Baker et al. 2003- The same 
procedure was used, but a fixed PSF derived from the 
unsaturated guide star in each SOFI science fields was 
provided to GALFIT instead of the model PSFs. Not all 
the galaxies detected in the NACO images are present 
in the SOFI data, but on the other hand the SOFI spa- 
tial coverage is ~ 5 times larger, so that in the end we 
detected a comparable number of objects: 181 in total, 
92 brighter than K s = 20. The obtained distribution is 
shown in Fig. ^] as a gray histogram. The distribution 
of R e shows similar behavior in both AO and seeing lim- 
ited data, suggesting that the observed distribution of R e 
in the SWAN fields is real and not due to selection ef- 
fects. Nevertheless, the NACO distribution shows an ex- 
cess of sources in the range R e = 0.4 — 1.2 with respect 



to SOFI, indicating that the better resolution of NACO 
is able to resolve many sources that in the SOFI data 
appear merely "compact" (e.g., a bright elliptical galaxy 
that NACO correctly recovers with n = 4 and R e — 1" is 
recovered by SOFI for a seeing of 0.7" as a compact source 
with R e = 0.1", while a disky galaxy with n = 1 and the 
same effective radius is correctly recovered by both) . 

Finally, the axis ratio distribution shows the co- 
sinusoidal distribution expected from random inclinations, 
except for a peak at bj a — due to very compact galaxies 
for which GALFIT could not calculate the axis ratio, and 
hence which were arbitrarily assigned to this bin. 

6.3. Completeness 

We have derived a first estimate of the completeness limit 
in our images using simulated galaxy profiles convolved 
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Fig. 11. The difference A<fi between the radial vector from 
the guide star to the source at offset (Ax, Ay) and the 
source's own major axis position angle 4> as measured with 
SExtractor (left panel) with no PSF correction and with 
GALF1T (right panel) using the model PSF. In each panel 
the distribution for objects with angular distance from the 
guide star < 17" and > 17" are shown (solid lines) as well 
as the total distribution of all sources (dashed line) . Error 
bars are yf~N for N the number of objects in each bin. 



with true NACO PSFs. This is a difficult task, as the 
completeness depends not only on the brightness and mor- 
phology of the source, but even on the distance from the 
guide star and the position on the frame. Therefore we de- 
rived a first estimate using only one reference PSF, a star 
at 21.2" from the guide star in SBSF 15. Fields with sim- 
ulated galaxies of fixed magnitude and different R e (from 
0.1" to 1.0") were created, and the sources were recov- 
ered using the same SExtractor parameters as used for 
the SWAN fields. We find that at K s = 20 all the simu- 
lated elliptical-like galaxies were detected, while the disk- 
like galaxies were only 100% complete for effective radii 
smaller than R e — 0.6". An average 50% completeness 
limit - i.e., the 50% limit for galaxies with R e = 0.3" 
- is K s ~ 21.5 (Kab ~ 23.3) for disk-like galaxies and 
K s ~ 23 (Kab ~ 24.8) for elliptical-like ones. A more 
detailed analysis of completeness will be presented in a 
forthcoming paper (Cresci et al. 2005). 

6.4. Confirming that we are accounting for 
anisopla na ticism 

The recovered position angle distribution can be used to 
indicate whether we are indeed correcting for the anisopla- 
natic PSF in an effective manner, and to evaluate possible 



biases in the results due to contributions by the wavefront 
correction that are still uncorrected by our PSF model. We 
have divided our sample in two subsets with separations 
AO < 17" and AO > 17". This separation was chosen em- 
pirically as the approximate radius beyond which the PSF 
distortion becomes noticeable (see Baker et al. 2005), but 
also corresponds to the average isoplanatic angle of the 
five fields in our analysis here. If no radial stretch due to 
PSF anisoplanaticism is present (or if this effect is fully 
taken into account by a PSF model) we would expect the 
distribution of A<f> = <f> — arctan(Ax/Ay), i.e., the differ- 
ence between the radial vector from the star to the source 
at offset (Ax, Ay) and the source's own major axis posi- 
tion angle <f>, to be uniform over the interval [—90°, +90°] 
for both AO subsamples. 

The distribution derived using SExtractor with no PSF 
anisoplanaticism correction (see Fig. ^] left panel) shows 
a clear peak at A(f) — for objects at large radial dis- 
tances from the guide star, i.e., these objects appear to be 
preferentially oriented towards the guide star. The right 
panel of Fig. II II shows the same distribution as measured 
by GALFIT using the model PSF. As one would expect, 
in this case the peak is greatly reduced due to the model 
PSF. The small remaining overcounts could easily be due 
to random chance rather than a systematic effect. Indeed, 
analysis of this effect in the simulated galaxies indicates 
that it is completely removed. 

7. Conclusions 

In this paper we have presented a new approach to account 
for the PSF variations across the field of view in deep 
AO images, developed to correct the anisoplanaticism in 
SWAN images obtained with NACO at the ESO VLT, but 
also more generally applicable to other wide-field adaptive 
optics data. The survey is intended to overcome the cur- 
rent shortage of extragalactic AO targets and to study 
faint and compact field galaxies with unprecedented reso- 
lution in the near-IR. The observing strategy fully exploits 
the present capabilities of AO instrumentation, so that our 
method for PSF modeling may be broadly useful for AO 
cosmology. 

We have described the PSF as the convolution between 
the on-axis PSF and a spatially varying kernel, an ellip- 
tical Gaussian elongated towards the AO guide star. We 
find that even adopting the most extreme case in which the 
the difference between the kernels in two different fields is 
given by a single parameter, namely the isoplanatic angle, 
the PSF can still be described with enough accuracy to 
extract reliable morphological parameters of field galax- 
ies. This approach is particularly convenient, as it uses 
only easily available data and makes no uncertain assump- 
tions about the stability of the isoplanatic angle during 
any given night. In addition, just a few point sources are 
required in each field to derive the isoplanatic angle and 
therefore compute the kernel. 

Our simulations demonstrate that the model is able 
to recover reliably morphological parameters in a typical 
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SWAN field with an integration time of one hour up to 
K s <~ 20.5 (Kab ~ 22.3), and that our very simple model 
is nearly as good as using the true PSF to account for 
anisoplanaticism. 

Finally, we have presented the first morphological re- 
sults for five SWAN fields, using the GALFIT package 
coupled with our model PSF. The recovered source pa- 
rameters confirm that we are indeed accounting for aniso- 
planaticism, and show the unique power of AO observa- 
tions to derive the details of morphology in faint galaxies. 
These results pave the way for the forthcoming analysis of 
all the obtained SWAN fields, which will combine for the 
first time the high angular resolution of a space-based sur- 
vey with the shallower depth and wider area of a ground- 
based survey. 
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